Junctions of three quantum wires and the dissipative Hofstadter model 
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We study a junction of three quantum wires enclosing a magnetic flux. This is the simplest problem 
of a quantum junction between Tomonaga-Luttinger liquids in which Fermi statistics enter in a non- 
trivial way. We present a direct connection between this problem and the dissipative Hofstadter 
problem, or quantum Brownian motion in two dimensions in a periodic potential and an external 
magnetic field, which in turn is connected to open string theory in a background electromagnetic 
field. We find non-trivial fixed points corresponding to a chiral conductance tensor leading to an 
asymmetric flow of the current. 
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Electric conduction in quantum wires is of much cur- 
rent interest. From a practical viewpoint, continuing 
advance in electronic technology is now reaching a level 
which requires understanding electric conduction in the 
quantum limit. From a basic science perspective, con- 
duction in low-dimensional systems exhibits many inter- 
esting properties due to strong correlation effects. A sem- 
inal example is discussed by Kane and Fisher (KF) 
who predicted that backscattering due to a single im- 
purity, however small, makes the conductance vanish at 
low temperature in the presence of repulsive interactions. 
This is a clear manifestation of the non-Fermi liquid be- 
havior of interacting electrons in one dimension, which 
can generally be regarded as a Tomonaga-Luttinger liq- 
uid (TLL). A TLL is nothing but a field theory of free 
bosons in 1 + 1 dimensions, and the effect of the impurity 
can be studied as a boundary interaction in the field the- 
ory; renormalization-group fixed points can be identified 
with conformally invariant boundary conditions. Study 
of boundary conditions in conformal field theory is an 
active research area in itself, with applications to string 
theory and statistical mechanics. 

In order to make any useful circuit with quantum wires, 
one needs a junction of three or more wires, which gener- 
alizes the KF problem. This turns out to be a surprisingly 
rich problem 0, 0, 0, IE 0] with a number of open 
questions, even for three wires. One of the novel aspects 
of the multi-lead problem is that, unlike the junction of 
two wires (KF problem), electron fermion statistics play 
a crucial role |2( ■ 

Here we study a junction of three quantum wires with 
an enclosed magnetic flux <f>, as shown in Fig. ^ Using a 
low energy limit Dirac fermion formulation, the physical 
junction becomes equivalent to electron hopping between 
wires, represented by the boundary term 
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H B = - ^[(re'Wtyt^ + h.c.) + rVj^-]. (1) 
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Here ipj is the electron annihilation operator for the wire 




FIG. 1: A quantum junction of three wires containing a mag 
netic flux d>. 



j at the junction (with y^ = ijj 3 ). T is the (real) hopping 
amplitude, r a reflection amplitude and <f> is essentially 
the dimensionless flux through the junction (for weak 
hopping), r has no effect on the conductance (in the zero 
frequency, zero temperature limit), with the exception of 
the non-interacting case, so we henceforth ignore it. For 
simplicity, we consider spinless fermions, while consider- 
ing arbitrary strength bulk interactions. In a TLL, the 
interaction strength is essentially contained in a single 
parameter <?, which determines various physical quanti- 
ties and exponents, g < 1 and g > 1 correspond, re- 
spectively, to repulsive and attractive interactions, while 
g = 1 is the Fermi liquid. 

In the presence of a magnetic flux <f> through the junc- 
tion, the combination of the quantum phases due to the 
flux and fermion statistics is mapped to a "magnetic 
field" for the free boson field at the boundary. Such a the- 
ory has been discussed from the viewpoint of open string 
theory 0. Moreover, the phase diagram in the pres- 
ence of both interactions and the "magnetic field" at the 
boundary was studied in a beautiful paper by Callan and 
Freed (CF) [j| . We will show that a low-energy renormal- 
ization group (RG) fixed point of our junction problem 
is given by the "magnetic" boundary condition of CF, 
which leads to a chiral conductance tensor correspond- 
ing to an asymmetric flow of the current. The "mag- 
netic" boundary condition exhibits a non-monotonic de- 
pendence on the TLL interaction parameter g, unlike in 
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other applications of TLL. In particular, the "magnetic" 
fixed point is stable for 1 < g < 3. This means that an 
arbitrarily small magnetic flux cj> leads to a substantial 
breaking of time reversal symmetry for this range of g. 

The electron annihilation operator in wire j is repre- 
sented as 



(2) 



where rjj is the so-called Klein factor satisfying {rjj, r]k} = 
25 jk, which is necessary to ensure the anticommutation 
relations of the fermion operators in different wires. The 
Klein factors may be represented by the Pauli matrices. 
In the low-energy limit, the (imaginary time) effective 
action for the three disconnected quantum wires is given 
by the three component free boson 



f 3 
S = I dr dx 



(3) 



defined over the half-line x > 0, with the junction at 
x = 0. Neumann boundary conditions dipj/dx = are 
satisfied. 

The boundary interactions lead to renormalization to 
an infrared fixed point with a different boundary con- 
dition. However, current conservation at the junction 
requires the "center of mass" field 



$0 = —7=(<Pl + ¥>2 + <f3) 
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(4) 



to always obey Neumann boundary conditions. Thus the 
remaining degrees of freedom at the boundary comprise 
a two component boson field $ = ($1,^2) defined by 
$1 = Oi - f 2 )/V2, $2 = (<fi + ¥2 - 2<p 3 )/V6. Using 
this new basis, the hopping term is written as 



f 3 

J 0=1 



JK a -<5> 



+ h.C, 



(5) 



where K x = (-1, V3)/2, K 2 = (-1, ->/3)/2, K 3 = (1,0). 
The scaling dimension of the hopping term in the dis- 
connected limit (r = 0) is calculated by standard meth- 
ods as 1/g, which is not affected by the Klein factors 
nor the magnetic flux. However, if one attempts a per- 
turbation theory in the hopping amplitude T, a phase 
factor appears in each term. This makes the problem 
rather different from the standard free boson field the- 
ory, and hence it is difficult to predict the nature of the 
strong hopping limit. In particular, for g > 1, when the 
disconnected wire fixed point is unstable, it is difficult to 
identify the infrared fixed point. 

Each order of the expansion of the partition function 
contains an integral over the correlation function 
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where 5k{L) = 1 if L = and 5k (L) = otherwise, Lj 
is one of the six vectors ±J?i.2,3, and ()o denotes the ex- 
pectation value with the Neumann boundary condition 
on <f>. Thus, identifying the vertex operator e lLj <s ' with 
the displacement by Lj in the two-dimensional plane, a 
non-vanishing contribution © corresponds to a closed 
loop on a triangular lattice spanned by K\ and K 2 . The 
perturbation series is given by the multi-point correla- 
tion function @ to any order, in principle. However, in 
the present problem, we pick up extra phase factors due 
to the Klein factors (and the flux <j>.) The phase factor 
for a given loop is given by the product of phase factors 
of elementary triangular loops. Referring to the two in- 
equivalent triangles in a triangular lattice as "up" and 
"down" , the phase factor for the counterclockwise loop 
on the up-triangle is determined to be 



(-m)e^ /3 (-iV2)e l<p/i (-ir,i)e 



i<f>/3 ( 
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i<t> 



(7) 



while that for the counterclockwise loop on the down- 
triangle is 



(tV3)e-^ /3 (iV2)e-^ /3 (m)< 



*^/3 — p i( T ' 



(8) 



Namely, loops on the triangular lattice pick up a phase 
as if there is a staggered magnetic flux of <f> and n — <f> in 
each elementary triangle. 

Fortunately, we shall see that the above perturbative 
expansion coincides with the "generalized Coulomb gas" 
studied by CF |9( in the context of the dissipative Hofs- 
tadter model, DHM, (quantum motion of a single particle 
under a magnetic field and a periodic potential, subject 
to dissipation.) The "free" action reads: 



S [X] 



1 



(9) 

where fx, v = 1,2 and e 12 — —e 2 \ = 1, en = e 2 2 = 0. 
a and /3 are related to the dissipation and the magnetic 
field, respectively. This determines the propagator: 



D^(r) = (JMt)X„(0)> = 



a 2 + f3 2 
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a 2 + p 



sgnre^. (10) 



We now introduce a potential term, which is somewhat 
different from the "rectangular" one in Ref. 0, as 



S V [X] = -Ve t5/3 / dr 



1^ 



C.C., 



(11) 



where V and 5 are chosen to be real. Expanding the par- 
tition function in powers of V , a nonvanishing contribu- 
tion is given by an integral over the correlation function 
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where Lj is again one of the six vectors ±^1^,3. 

The perturbation series in the DHM can be made to 
match that in the quantum wire model, including the 
Klein factors and flux, for appropriate choice of a and 
[3. The absolute value of the correlation functions agree 
exactly if 
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The phase factor for any loop in the DHM is determined 
by that for two elementary triangles: 
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±5, (14) 



where + and — signs apply to the up and down triangles, 
just as in the quantum wire problem, 10 and JHJ. Thus, 
if one chooses 



V3n- 



P 



= (2n- 1)tt 



(15) 



with an integer n, there is an appropriate 5 which makes 
the two expansions based on ©-© and (|12|l coincide 
exactly including the phases. Because the phases are only 
defined modulo 2tt, there is actually an infinite number of 
different choices of a and (3 labeled by the integer n. Each 
choice defines a quite different theory with respect to the 
dynamical variable X, but gives the identical generalized 
Coulomb gas. 

Now let us consider the case of g > 1. Then, as we have 
already discussed, the electron hopping is a relevant per- 
turbation at the disconnected limit. In the DHM, for any 
choice of n, V is a relevant perturbation for g > 1. We 
would like to find the infrared stable fixed point reached 
in the low energy limit. A simple guess is that it occurs 
at V — > 00. The stability of the V — > 00 fixed point can 
be determined using the instanton method [|J. Namely, 
in the strong potential limit, the X field is pinned at one 
of the minima of the potential (f 1 If) . The leading pertur- 
bation to this limit is given by a tunneling between the 
neighboring minima, represented by an "instanton." A 
calculation following Ref. 9] gives the scaling dimension 
of such perturbation to be a times the distance between 
the nearest minima squared where a is determined by 
p3 [l -i|15 [l . For a generic value of (f>, minima of the poten- 
tial Hll|) form a triangular lattice with lattice constant 
2/y/3. Thus the dimension of the most relevant operator 
at this V — * 00 fixed point is: 
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3+ (2n- 1)V 



(16) 



This implies that the V — > 00 fixed point is unstable for 
any value of g for all choices of n except for n = or 1 



where it is stable for 1 < g < 3. In this range of g we 
expect that the infrared limit of the junction model cor- 
responds to the V — > 00 limit of the DHM for n = or 1. 
The infrared fixed point for other values of n must corre- 
spond to nontrivial intermediate V fixed points, for which 
we do not know explicit solutions. The infrared fixed 
point for each value of n should give the same physical 
behavior for the junction model. In principle we could 
study this using any choice of n. However, clearly the 
n = or 1 choices are preferred because only in those 
cases can we explicitly analyze the fixed points. Let us 
denote the RG fixed points corresponding to the strong 
potential limit in these cases as \- an d X+j respectively 
for n — and 1. Since x± have the same stability, at 
this stage we cannot determine which of these two fixed 
points is realized in the infrared limit for a generic value 
of <f>. 

On the other hand, in the special cases <j> — ±7r/2 
which maximally break the time reversal symmetry, the 
choice of the fixed point is unique. For <p — tt/2, the 
choice n = 1 (P positive) gives 5 = 0, and the poten- 
tial minima forms a triangular lattice as usual. However 
the other choice n = (P negative) gives 6 = ir, for 
which the potential minima form a honeycomb lattice 
with minimum distance 2/3, making the strong poten- 
tial limit unstable for all g. Similarly, for <f> = —ir/2, 
n = (P negative) is the unique choice to give a stable 
fixed point. This suggests that these fixed points \± re_ 
fleet the breaking of time reversal symmetry due to the 
magnetic flux cf>. Indeed, the conductance at these fixed 
points exhibits a chiral behavior breaking the time re- 
versal invariance, as we show below. On the other hand, 
changing the flux corresponds to an irrelevant perturba- 
tion at these fixed points. Thus we conjecture that the 
flow from small V goes to the V — > 00 fixed point in the 
n = 1 representation, but not in the n — representa- 
tion, where we expect the flow to be to a non-trivial fixed 
point, for < <j> < it. Conversely, for — it < <j> < 0, we 
conjecture that the infrared fixed point is given by the 
V — > 00 one in the n = representation, not n = 1. 

The conductance at the junction is expressed as a ten- 
sor 



(17) 



where Ij is the total current flowing into the junction 
from wire j and Vk is the voltage applied to wire k. From 
current conservation Gjk = 0. It is also clear that the 
currents depend only on voltage differences. According 
to the Kubo formula, the conductance tensor is given by 



G 



jk 



lim 



1 



d 2 lnZ 



w^o+ uj dAj(— u>)dAh(u>) ' 



(18) 



where u> is the Matsubara frequency, and Aj and Vj are 
related in real time t by dAj/dt = —Vj. Assuming that 
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the voltage drops only across the junction, its effect can 
be represented as an extra time-dependent phase factor 
(vector potential) e ltyA i~~ Ak ^ in the electron hopping term 
from wire k to wire j. It can be conveniently included in 
the CF representation as 



S V [X] = -Ve l5/3 



dTj2e iKr(X+d) 

3=1 



C.C. 



(19) 



where a = J2j k i Kj€jki(Ak — Ai)/3. Thus the conduc- 
tance is related to 



d 2 hxZ 
dAjdAk 



£ e 3m e kn K X ^^ (20) 



where e 12 = e 2 3 = £31 = 1, £21 = £32 = £13 



-1 and 



ejj = 0, /i, v = 1,2 refers to the components of two- 
dimensional vectors a and Kj. In fact, 

d 2 InZ 



da>*{-bJi)da v (w) 



(pJ<J)p v (-u)) 



(21) 



is the correlation function of the generalized Coulomb gas 
density as defined in Ref . 9] , which we use to determine 
the conductance tensor. In the strong potential limit 
of the CF representation, the X field is pinned and the 
correlation function is given exactly by 



\XV ■ 



(22) 



From (|18fl - l|22|). we obtain the conductance tensor 
Gf k = G[(3S jk -l)±ge jk ]/2, 



(23) 



with G = (e 2 / h)4g / (3 + g 2 ) and ± corresponding to the 
fixed points x±- The term proportional to ey/., repre- 
sents an asymmetric conduction which is allowed by the 
Z3 symmetry of the Hamiltonian once time-reversal sym- 
metry is broken by the flux. 

For example, when a voltage V\ > is applied to wire 
1 while V2 = V3 = 0, the current flowing from wire 1 to 
the junction is I± — GV±, which is independent of the 
asymmetry. On the other hand, the currents flowing out 
from the junction to wires 2 and 3 are given by — 12 = 
G(l ± g)V x /2 and -I 3 = G(l T g)Vi/2. In fact, for g > 
1, the current for fixed point x+ (X— ) rather flows in 
towards the junction from wire 3 (wire 2), opposite to 
the voltage drop, when the voltage V% > is applied 
only on wire 1! 

The asymmetric conduction is forbidden under time 
reversal but is possible in the presence of magnetic flux 
4>. In the limit of g — » 1, which corresponds to non- 
interacting electrons, we find for example Gn = — G\2 = 
lj G13 = for x+- This asymmetric conduction of 
non-interacting electrons can be simply understood as 



a boundary condition •tp™ 
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At g = 1, there is 



in fact a continuous family of boundary conditions char- 
acterized by a 3 x 3 S-matrix S as ip° ut = J2 k Sjk^\ n ■ 



Our x± correspond to special cases of these free fermion 
boundary conditions in the g — > 1 limit. 

A remarkable aspect of the fixed points x± is that the 
conductance l|23|l . as well as the scaling dimension (TSJ, 
exhibits non-monotonic dependence on the interaction 
parameter g, unlike in other known applications of TLL. 
Making the electron interaction more attractive can de- 
crease the conductance! 

To summarize, we studied a junction of three quan- 
tum wires of interacting spinless electrons with a mag- 
netic flux. The fermionic statistics, together with the 
magnetic flux, bring nontrivial phases into the problem, 
making it different from the standard boundary problem 
of free boson field theory. It is shown to be equivalent 
to a certain generalized Coulomb gas introduced by CF. 
Using the mapping, we have shown that the chiral fixed 
points x± exhibit an asymmetric conduction and non- 
monotonic dependence on the interaction parameter g. 
The fixed points x± are stable for 1 < g < 3, and we ex- 
pect the junction with <fi = tt/2 (<f) — —tt/2) to be renor- 
malized to the fixed point x+ (X-)i since it is the unique 
stable fixed point obtained with the CF representation 
at V — > 00. We conjecture that the system renormalizes 
to these fixed points (for 1 < g < 3) for all values of (j> 
except those which respect time-reversal: <j>/n integer. 

We note in passing that our results also apply to an 
equivalent problem. If we interpret $ as the co-ordinate 
of a neutral spin-1/2 particle, 10 and JSJ describe its 
propagation in a periodic textured Zeeman field with 
components B a (Q) — 4Tsm(K a ■ $ + ^/3), and Hamilto- 
nian Hz — —!?($) • S, subject to dissipation. 
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